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ABSTRACT 

Using cosmological A^-body simulations of dark matter haloes, we study the effects 
of non-sphericity, substructure and streaming motions in reproducing the structure and 
internal kinematics of clusters of galaxies from kinematical measurements. Fitting an 
NFW model to the 3D density profile, we determine the virial mass, concentration 
parameter and velocity anisotropy of the haloes, and then calculate the profiles of 
projected velocity moments, as they would be measured by a distant observer. Using 
these mock data, we apply a Jeans analysis for spherical objects to reproduce the linc- 
of-sight velocity dispersion and kurtosis and fit the three parameters. We find that the 
line-of-sight velocity dispersion and kurtosis profiles of a given halo can vary consider- 
ably with the angle of view of the observer. We show that the virial mass, concentration 
parameter and velocity anisotropy of the haloes can be reproduced satisfactorily in- 
dependently of the halo shape, although the virial mass tends to be underestimated, 
the concentration parameter overestimated, while the recovered anisotropy is typically 
more tangential than the true one. The mass, concentration and velocity anisotropy of 
haloes are recovered with better precision when their mean velocity profiles are near 
zero. 

Key words: methods: iV-body simulations - methods: analytical - galaxies: clusters: 
general - galaxies: kinematics and dynamics - cosmology: dark matter 



1 INTRODUCTION 

There is a long tradition of determining the internal kine- 
matical properties of bound systems based on the Jeans 
equations, which are velocity moments of the collisionless 
Boltzmann equation, and which link the 3D velocity mo- 
ments (e.g. velocity dispersion and kurtosis) to the poten- 
tial gradient. The Jeans equations are thus used to model 
the mas s and velocity distribut ion inside elliptical galax- 
ies (e.g. |Bjrmey & Mamon 1982), clusters of galaxies (e.g. 
iKent fc GunrJll982l:lMerritdll987|h as well as globular clus- 
ters (e.g. iMerrifield fc KentHl990ft . The standard Jeans ap- 
proach for the determination of the mass and velocity dis- 
tribution assumes equilibrium and sphericity of the system. 
However, even in the inner parts of an object, there is non- 
virialized matter, for example matter falling into a cluster for 
the first time or smaller clumps in the process of relaxation. 
Such non-virialized matter could produce misleading results 
when a cluster is studied through a Jeans analysis. Moreover, 



clusters of galaxies are not observed to be s pherically sym- 
metric feinggelilll983lWang fc Ufmerlll997lh nor are simu- 
lated structures of dark matter particle s with the masses o f 
clusters of galaxies dCole fc LacevHl99&|jiiig fc Sutdl2002lh 

One way to avoid non-virialized matter within a cluster 
is to restrict the Jeans analysis to the population of ellip- 
tical galaxies, which is thought to be dynamically relaxed 
jTullv fc Shavalll984t iBiviano et aHl2003l : lLokas fc MamorJ 
120031) . The question of how the existing substructure and 
non-sphericity may affect the results can only be fully ad- 
dressed by cosmological iV-body simulations including real- 
istic galaxy formation, where all 3-dimensional information 
would be available. 

The effect of incomplete virialization of structures of 
dark matter particles seen in cosmological iV-body sim- 
ulations on the estimates of the mass of a single clus- 
ter through the Jeans equation has been addressed by 
iTormen. Bouchet fc White! dl997l) . They showed that even 
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for significantly perturbed haloes, the mass M(r) at dis- 
tances larger than 2% of the virial radius inferred by the 
proper Jeans analysis is within 30% (r.m.s.) of the true mass 
and departs from it by less than 20% (r.m.s.) for average or 
relaxed haloes. 

In this work, we use cosmological iV-body simulations 
and analytical modelling to study the effect of departures 
from equilibrium and non-sphericity of dark matter haloes 
on the inferred properties of the halo and velocity distri- 
bution of its particles. We measure the mass and velocity 
distribution in the haloes and calculate the projected veloc- 
ity moments as an observer would do. We then perform a 
kinematic analysis based on the Jeans equations to check to 
what extent we can reproduce the properties of the haloes 
from the second-order (line-of-sight velocity dispersion) and 
fourth-order (line-of-sight kurtosis) velocity moments. 

The paper is organized as follows. In Section 2, we de- 
scribe the iV-body simulations used and calculate the prop- 
erties of the haloes chosen for analysis. In Section 3, we esti- 
mate projected velocity moments of the haloes. Section 4 is 
devoted to analytical modelling of those moments based on 
Jeans formalism and testing its reliability in reproducing the 
properties of the haloes. The discussion follows in Section 5. 



2 THE SIMULATED DARK MATTER HALOES 

We have used the iV-body simulations carried out by 
lHatton et all (l2003t) with their GallCS hybrid iV-body/semi- 
analytic model of hierarchical galaxy formation. The de - 
scription of this model can be found in lHatton etafl J20031) . 
The iV-body simulation contains 256 3 particles of mass 
8.272 x 10 9 M© in a box of size 150 Mpc and its softening 
length amounts to spatial resolution of 29 kpc. The simula- 
tion was run for a flat universe with cosmological parame- 
ters D.0 = 0.333, Q A = 0.667, H = 66.7 km s _1 Mpc -1 , and 
as = 0.88. Once the simulation is run, haloes of dark mat- 
ter are detected w ith a 'friends-of-friends' (FOF) algorithm 
iDavis et al.lll985T) using a variable linking length such that 
the minimum mass of the FOF groups is 1.65 x 10 11 M Q (20 
particles) at any time step. With this method, over 2 x 10 4 
haloes are detected at the final timestep, corresponding to 
the present-day (z — 0) Universe. 

We restrict our analysis to dark matter particles. Al- 
though the GallCS simulations we use include galaxy for- 
mation and evolution, this part of the simulations is based 
on a semi-analytical approach, which for our purposes is 
not yet satisfactory. For example, the isotropic velocity dis- 
tribution of GallCS galaxies is imposed and is not a result 
of virialization. We therefore conclude that the galaxies in 
these simulations are not reliable tracers of the overall dy- 
namical properties of the haloes and cannot be used to infer 
the density and velocity distributions from the 'observed' 
velocity moments. 

For our analysis, we have chosen the ten most mas- 
sive haloes formed in the simulation box (labelled here- 
after in order of decreasing virial mass as halo 1, halo 2 
and so forth). We estimated their virial radii, rioo, as the 
distances from the centre where the mean density is 100 
times the present critical density (i n agreement with the so- 
called spherical collapse model, see lKitavama fc Suto Il996t 
lLokas fc Hoffman! l200ll) . The centres of the haloes are de- 



Table 1. Three-dimensional properties of the simulated haloes 
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2.3 


1.3 


1 


16.7 


3.2 


5.6 


0.40 ± 0.06 


0.005 


2 


1.9 


1.1 


1 


8.5 


2.5 


5.0 


0.13 ± 0.03 


0.103 


3 


2.8 


2.0 


1 


8.1 


2.5 


4.5 


0.24 ±0.05 


0.134 


1 


3.3 


1.3 


1 


7.0 


2.4 


7.1 


0.14 ±0.06 


-0.260 


5 


2.0 


2.0 


1 


4.5 


2.1 


7.9 


0.15 ± 0.07 


0.065 


6 


2.6 


1.4 


1 


1.3 


2.0 


5.0 


0.36 ± 0.05 


0.028 


7 


1.8 


1.4 


1 


4.1 


2.0 


8.9 


-0.03 ± 0.05 


-0.079 


8 


2.4 


1.6 


1 


3.9 


2.0 


10.0 


0.44 ± 0.03 


-0.018 


9 


2.2 


1.5 


1 


3.2 


1.8 


7.1 


0.08 ± 0.12 


-0.012 


10 


1.8 


1.3 


1 


2.9 


1.8 


10.0 


-0.23 ± 0.08 


-0.168 



termined as the local density maxima which turn out to be 
slightly different than the centres of mass found using the 
FOF algorithm. Within the virial radius, the haloes have 
2 x 10 s particles for the most massive halo and 4 x 10 4 par- 
ticles for the least massive one of the haloes we have chosen. 
The masses of the haloes, M100 = M(rioo), are listed in 
Table □ 

All haloes have similar 3D radial phase space distri- 
butions. An example of such a distribution for the most 
massive halo is presented in Figure Q The Figure shows 
the radial velocity (with respect to the centre of the halo) 
of the particles in u nits of the c ircular velocity at rioo 
(wioo = fc(^ioo) = \f GMioo/rioo) as a function of radial 
distance measured in units of rioo- We have divided the 
phase space into different regions corresponding to differ- 
ent dynamical states that we expect to find in a typical 
dark matter halo. Inside the virial radius, we can identify 
three regions. The region with low to moderate absolute 
velocities is likely to be populated by virialized particles, 
although there are indications of a group of particles at 
r = 0.3 rioo slowly moving into the halo core. On the out- 
skirts of the velocity distribution we can find particles whose 
dynamical state is not clear, as they could be either high- 
velocity outliers of the virialized component or else infalling 
towards the cluster core or already in a rebound regime after 
a passage through the centre. Beyond the virial radius, we 
can also find particles in a rebound trajectory and particles 
on the infalling branch (which includes particles expand- 
ing away from the cluster beyond the turnaround radius at 
2.5-4 virial radii). In Figure these five different subsam- 
ples are denoted by virialized, rebound, inf all, vir-inf ? 
and vir-reb?, respectively, where the question marks indi- 
cate the uncertainty in the actual dynamical state of these 
subsamples. 

It is important to note that when studying projected 
quantities we need to consider all five subsamples, as well as 
the whole sample (which will be marked hereafter by 'all'), 
because, as seen projected on the sky, there are particles be- 
longing to each subsample that fall inside the 'virial' cylinder 
whose projected radius is the virial radius rioo, but there is 
no way of determining which of these particles are actu- 
ally within the 'virial' sphere of radius rioo (inscribed in the 
virial cylinder). 

Figure|2|shows the density distribution in the most mas- 
sive halo 1 as a function of radial distance in units of rioo. 
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Figure 1. 3D radial phase space diagram (radial velocity, Hub- 
ble flow included, versus radial distance, both with respect to 
the cluster centre and normalized to virial quantities fioo an d 
rioo, respectively) for a random set of 1% of particles of the most 
massive of the simulated haloes (halo 1). The oblique lines sep- 
arating the dynamical regimes were drawn 'by eye' and corre- 
spond to v r /vioo = 1-8 — 1.06r/rioo (bottom) and Vr/vioo = 
-1.8 + 1.06 r/noo (top). 
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Figure 2. Density profile of halo 1. The measurements were done 
in radial bins of equal logarithmic length and errors were esti- 
mated as Poisson fluctuations. The error bar in the second point 
shows the maximum error assigned to measurements. The curve 
is the best-fitting NFW profile with concentration c = 5.6. 



We found that the measure d density profile is well approxi - 
mated by the NFW formula (iNavarro. Frenk fc Whitell997l) 



Q(s) = A e c 2 ff(c) 
Q c ,o 3s(l + cs) 2 



(1) 



where s = r/r^o, Q c ,o is the present critical density, A c = 
100, c is the concentration parameter and g(c) — [ln(l + c) — 
c/(l + c)] _1 . The statistical errors are much smaller than the 
departures due to substructure in the halo, but the overall 
fit is satisfactory. We find that the best-fitting concentration 
parameter for halo 1 is c = 5.6. The virial mass of halo 1 is 
Mioo = 1.67 x 10 15 M Q so the concentration we estimated 
is consistent with the d ependence of c on ma ss inferred from 
iV-body simulations bv lBullock et all j2001), also run with 
a ACDM cosmology. Similar results are obtained for other 
haloes. The fitted parameters are summarized in Tabic Q 

In Figure El we present radial profiles of the mean ra- 
dial velocity in units of the circular velocity at rioo and the 
anisotropy parameter 



/3 : 



*f(r) 
cr?(r) 



(2) 



for the ten haloes, where ag and o> are the velocity dis- 
persions (with respect to the mean velocities) discussed in 
detail in the next Sections. We show measurements for all 
particles inside the sphere of radius rioo- The curves rep- 
resent the anisotropy or mean radial velocity of dark mat- 
ter particles enclosed in shells of thickness O.lrioo centered 
at 0.05rioo, 0.15rioo etc. As we can see, the orbits of most 
of the haloes are mildly radial, with positive mean f3. We 
have calculated the unweighted mean anisotropy inside the 
virial radius. These values are listed in Table Q where the 
uncertainties are the dispersions of the values about the un- 
weighted mean. The anisotropy does not vary strongly with 
radius so in the modelling which follows we will assume it to 
be constant and equal to the mean value inside the virial ra- 
dius. However, it is interesting to note that the variation of 
(3 with distance from the centre is very different for each of 
the analyzed haloes. As for the mean radial velocity, we can 
see from Figure [3] that it is consistent on average with zero 
inside rioo for most haloes, albeit with non-negligible radial 
variations caused by internal streaming motions, including 
some abrupt variations suggesting the presence of smaller 
haloes (akin to groups of galaxies) falling in or bouncing 
out of the halo. In Table we list this quantity in units of 
wioo- The negative sign indicates infall motion towards the 
centre of the halo. 

Since we wish to study the effect of non-sphericity of 
haloes by choosing different directions of observation, we 
have determined the principal axes of our haloes from their 
moments of inertia, using particles at radial distances up to 
rioo - The ratios of the eigenvalues of the inertia tensor, listed 
in Table show significant departures from sphericity. 



3 VELOCITY MOMENTS OF DARK MATTER 
PARTICLES 

We now study the kinematical properties of haloes as they 
would be seen by a distant observer. The quantities dis- 
cussed are all projected along the line of sight. With the 
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Figure 3. The radial profile of the anisotropy parameter /3 (top 
panel) and the mean radial velocity in units of the circular velocity 
at rioo (bottom panel) for all particles inside rioo for the 10 haloes 
of Table 

Jeans formalism, we can only model the projected veloc- 
ity moments of objects in equilibrium, and thus we cannot 
make any prediction about the velocity moments of the re- 
gions of phase space called 'rebound', 'infall', 'vir-reb?' and 
'vir-inf?' in Figure Q Moreover, not all particles inside rioo 
are virialized, as they could be falling into the core for the 
first time or be on a rebound orbit not yet in equilibrium. 
However, the division of the phase space into the 5 regions 
(Fig.0 was made by eye, with no exact determination of the 
dynamical state of the particles, and it was only based on 
the radial component of the particle velocities. This could 
produce misleading results in the analysis of the velocity 
moments, which involve all components of the velocity. For 
these reasons, in what follows, we will restrict ourselves to 
the study of the velocity moments of all particles inside the 
virial sphere of radius rioo and all particles within the virial 



cylinder of projected radius smaller than rioo- The latter 
would be the ones used in the Jeans formalism by an ob- 
server unable to distinguish which particles actually lie in 
the virial sphere of radius rioo- 

We mimic the observations as follows. For each halo, we 
place an observer at either 0°, 45° or 90° with respect to the 
major axis (see Tabled so that the three chosen directions 
are in the same plane defined by the major axis and the sum 
of the other two axes. The choice of axes is dictated by the 
non-sphericity of the set of particles within the virial sphere, 
but also of the filaments of groups and other matter falling 
into the virial sphere out to scales ~ 10 rioo- It turns out 
that the principal axes on scales of ~ 10 rioo are very similar 
to those computed for r < rioo, so we restrict ourselves to 
the principal axes obtained from the tensor of inertia of the 
particles within the virial sphere. We then project all the 
particle velocities along the line of sight and the distances 
on the surface of the sky. In Figure 0] we show line-of-sight 
velocities of dark matter particles in halo 1 projected along 
the major axis as a function of projected distance from the 
centre in units of rioo- 

Observers remove a fraction of the interlopers of a clus- 
ter by excluding the high-velocity outliers. Here, we remove 
the high- velocit y outliers from our mo ck sa mples in a si milar 
way a s done bv lKent fc Gunnl l)l982f) and lLokas fc Mamonl 
( 2003) for clusters of galaxies. The velocity cuts are shown as 
solid curves in Figure 0] Less than 2% of all particles within 
the virial cylinder were removed in this fashion. Although 
the main body of the halo is quite well defined in velocity 
space, with much more particles than there are galaxies in a 
cluster, the gaps between t he halo and the backgro und are 
not as visible as found bv lLokas fc Mamonl j200ot) for the 
Coma cluster. After applying this selection procedure, the 
fraction of particles lying inside the cylinder of projected ra- 
dius rioo that are actually outside the sphere of radius rioo 
is between 6% and 35%, with a mean of 15%. 

We divide the projected radius in ten bins and calcu- 
late, in each bin of projected radius, the mean, dispersion, 
skewness, and kurtosis of the line-of-sight velocities, Vi, ac- 
cording to the following formulae 
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where N represents the number of particles per bin. 

To see the effects of non-sphericity, we show the results 
for haloes 1, 4 and 5, which have three different shapes: 
halo 5 is oblate, halo 4 is roughly prolate and halo 1 is tri- 
axial. The results for halo 1 are shown in Figures |3J7| for 
observers situated at 0°, 45° and 90° with respect to the 
major axis. The errors were estimated using bootstraps, but 
since the number of particles in each bin is very large (of 
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Figure 4. Linc-of-sight velocities of a random set of 2% of the 
dark matter particles in halo 1 projected along the major axis as 
a function of projected distance from the centre in units of rioo- 
The curves indicate the cuts used to distinguish between particles 
belonging to the halo from the background. The vertical clumps 
indicate groups of particles (nearly all within the virial sphere). 
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Figure 5. Projected velocity moments of the dark matter parti- 
cles in halo 1 measured at 0° to the principal axis. The upper left 
panel shows the mean line-of-sight velocity with respect to the 
velocity of the centre of the halo in units of the circular veloc- 
ity at rioO' The upper right panel gives the line-of-sight velocity 
dispersion in the same units. The two lower panels give the skew- 
ness (left) and kurtosis (right). In each panel the solid line shows 
results for particles lying inside the sphere of radius riod while 
the dashed line is for all particles. The dotted curve shows the fits 
obtained from the Jeans equations. 
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Figure 6. Same as Figure|2] but measured at 45° to the principal 
axis. 



Halo 1 : 90 degrees to principal axis 
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Figure 7. Same as Figure|S] but measured at 90° to the principal 
axis. 



the order of 10 4 ), the errors are small and do not account 
for the variability of the profiles, which is mainly due to 
substructure. In each panel, the dashed line shows results 
for all particles that in projection end up with a projected 
radius smaller than the virial radius, whether or not they 
are actually within the virial radius in 3D space. The solid 
line shows results for the particles really lying inside rioo- 
In the right-hand panels presenting even moments (velocity 
dispersion and kurtosis) in Figures EH3 we also show dotted 
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Figure 8. Same as Figures l5l7l but for halo 4. Line-of-sight skew- 
ness is not plotted. No theoretical predictions are shown. 



lines resulting from the fitting procedure based on the Jeans 
formalism presented in Section 4. 

The variation of the mean velocity with projected ra- 
dius provides an indication of the amount of substructure 
present in the halo. For observers at 0° and 45° with re- 
spect to the major axis of halo 1, the mean velocity with 
respect to the centre of the halo is approximately zero for 
every radial bin, indicating, if not the lack of substructure, 
at least the compensation of effects of different substruc- 
tures. This is true for the two subsets of particles studied, 
especially for particles actually within the virial sphere. For 
an observer at 90° with respect to the major axis of halo 1, 
we find a departure of the mean velocity from the velocity 
of the centre of the halo, which indicates the presence of 
substructure. As could be expected, this radial variation of 
the mean velocity is more pronounced for the particles we 
find in projection inside rioo that are not necessarily within 
the virial sphere of radius rioo- 

Contrary to the case of elliptical galaxies, where veloc- 
ity moments are measured from spectra obtained in slits, 
e.g. along the major axis, and thus do not involve averaging 
in shells of similar projected radius, here and in the analysis 
of galaxy motions in clusters the mean projected velocity is 
less affected by the presence of global motions like rotation 
or infall. Therefore, the estimation of the bulk velocities re- 
quires the 3D information. We find that the mean velocities 
in the radial and tangential directions inside the virial radius 
are typically of the order of few percent of uioo, except for 
halo 4 and halo 10 where significant net radial infall seems 
to be present (see Table 0. 

Figures |H| and [5] illustrate how the velocity moments 
depend on the shape of the given halo. To save space, we 
only show the even moments that we will model in the next 
Section and the mean line-of-sight velocity with respect to 
the velocity of the centre to test the relaxation of the haloes. 
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Figure 9. Same as Figure 181 but for halo 5. 



For the sake of clarity of the plots, we do not show the fits 
based on the Jeans equations, although we will discuss them 
in the following section. 

Even with significant noise from substructure, such as 
the trend for high mean velocity at R ~ 0.5 rioo for halo 
5 (Fig. 9) probably arising from a group bouncing out of 
the halo, we can see some common trends in Figures |SH5] 
First, the line-of-sight velocity dispersion and kurtosis pro- 
files can differ substantially for a given halo seen along three 
different axes (especially in halo 4). Moreover, the velocity 
moments are not strongly affected by the presence of parti- 
cles outside the virial radius, i.e. the solid and dashed lines 
typically do not differ significantly. This is understandable 
as the virial sample is a subsample of the 'all' sample and the 
mean difference in the number of particles is less than 20% 
as previously noted. The discrepancies between the velocity 
moments are more important for larger projected distances 
from the centre of a halo, as the surface density of the 3D 
haloes decreases faster than that of the surrounding mate- 
rial. 



4 MODELLING OF THE VELOCITY 
MOMENTS 

We now briefly describe the Jeans formalism for mod- 
elling the projected velocity moments of virialized objects 
and apply it to recover the properties of the haloes. The 
detailed description of t he calculations involved can be 
found in lLokasI i 2002l) andlLokas fc Mamonl 2003J) (see also 
iMerrifield fc Kentlll99"ol : Ivan der Marel et alJl200dl . 

Our purpose here is to reproduce the projected velocity 
moments discussed in the previous Section by solving the 
Jeans equations for the second and fourth velocity moments 
and adjusting the parameters describing the mass and ve- 
locity distribution in the haloes. We will then verify whether 
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the fitted parameters match the real properties of the haloes. 
The Jeans analysis assumes that the system is spherically 
symmetric and in equilibrium, that there are no net stream- 
ing motions (no infall and no rotation) so that the odd veloc- 
ity moments vanish. As we have seen in the previous Section, 
none of these is exactly the case for dark matter haloes. We 
want to check to what extent violating these assumptions 
affects the recovered properties of the haloes. 

One difference with respect to the analysis of galaxies is 
that the dark matter particles are very numerous in our sim- 
ulations (of order 10 5 per halo) while the number of galaxies 
in a cluster usually does not exceed a thousand. This is the 
reason why the errors due to substructure will be more sig- 
nificant here than sampling errors which are the dominant 
ones in measured velocity moments of galaxies. 

The second order velocity moments are v 2 and v 2 = i> 2 
and we will denote them hereafter by cr 2 and af. respectively. 
The y can be calculated from the lowest order Jeans equation 
(e.g. iBinnev fc Mamonlll982ft 



d . „. , 2/3 2 d$ 

— (ya r ) H va r = -v— , 

dr r dr 



(7) 



where v is the 3D density distribution of the tracer popula- 
tion and $ is the gravitational potential. Since in our case 
dark matter particles trace their own gravitational poten- 
tial, we have u(r) = g(r). We assume that the dark matter 
distribution is given by the NFW profile Q characterized 
by its virial mass Mioo and concentration c. We solve equa- 
tion Q assuming the anisotropy parameter of equation 
to be constant with — oo < f3 ^ 1. This model covers all in- 
teresting possibilities from radial orbits (/3 = 1) to isotropy 
(/3 — 0) and circular orbits (/3 — > — oo). 

The solution of the lowest order Jeans equation with 
the boundary condition a r — > at r — > oo for (5 = const is 
(e.g. lLokas fc Mamonll2003ft 



v<j r (j3 = const) = r 



-2/3 



2/3 d$ 

r v dr 

dr 



(8) 



As discussed in the previous Section, the quantity an ob- 
server would measure is the line-of-sight velocity disper- 
sion obtained from the 3D velocity dispersion by integrating 
along the line of sight iBinnev fc Mamonlll982l) 



I(R) 



\[r 



R 2 



dr 



(9) 



where I(R) is the surface distribution of the tracer and R 
is the projected radius. In our ca se I(R) is given by th e 
projection of the NFW profile (see lLokas fc Mamonll20oil) . 
Introducing equation Jj5J into equation @ and inverting the 
order of integration, the calculation of ai oa can be reduced to 
one-dimensional numerical integration of a formula involving 
special functions for arbitrary f3 = const. 

It has been established that systems with different den- 
sities and veloci ty anisotropies can produce identical tr i OB (-R) 
profiles (see e.g. lMerrifield fc Kentll99dlMerrittll987ft . This 
degeneracy can be partially lifted through the modelling of 
the fourth-order moment. With f3 = const, the solution of 
the Jeans equation for the fourth-order moment 



2/3.. -r 



I 41 i "1" 4 I o 2<1$ 
-(VVr) H VVr + &VO r 







dr r dr 

is (see lLokadl2002t lLokas fc Mamonll2003l) 



(10) 



fVr(/3 — const) = 3r 



-2/3 



2/3 2, 

r va r (r)—— dr 
dr 



(11) 



By projection, we obtain the line-of-sight fourth moment 
2 



< S (R) 



I(R) 



sfr 



R 2 



g(r,R,p) dr 



(12) 



where g(r,R,/3) = 1 - 2f3R 2 /r 2 + /3(1 + /3)i? 4 /(2r 4 ). 

Introducing equations ||SJ and (1111) into 112H and invert- 
ing the order of integration, the calculation can be reduced 
to a double integral. In the following, we use the fourth 
projected moment scaled with <rf OB in the form of projected 
kurtosis 



s(-R) 



,{R) 



(13) 



where the value of ki os (R) = 3 valid for a Gaussian distri- 
bution has been subtracted. 

We can now calculate the predictions of equations 
@ and 11211 for a given mass distribution and velocity 
anisotropy j3. As already discussed, we assume that the mass 
is given by the NFW distribution. For halo 1 we found the 
concentration of c = 5.6. For this value and different we 
obtain the profiles of <ti os (_R) and K\ os (R) shown in Figurc tTUl 
The values of a\ os are expressed in units of wioo and distances 
are in units of rioo- With these scalings, the predictions do 
not explicitly depend on the mass of the halo. The lines show 
results for different values of f3 = —0.3, 0, 0.3, 0.6 and 0.9, as 
indicated. 

We see that for increasingly rad ial orbits (in creasing f3) , 
the profile of <ti os turns steeper (e.g. lTonrvligSSh . Moreover, 
the kurtosis profile becomes more convex for increasingly 
radial orbits as opposed to the concave shapes in the case 
of isotropic and circular orbits. Since our measured kurtosis 
profiles in the previous Section have a concave shape and are 
slightly negative we do not expect the orbits to depart signif- 
icantly from isotropic, which is consistent with the measured 
anisotropy parameter (see Table 0. 

Mimicking the procedure used by observers to infer the 
mass and anisotropy profiles of galaxies and clusters, we fit 
the measured profiles of ai os and ki os of the ten haloes by 
solving equations JUJ and 11211 and adjusting the parameters 
Mioo, c and j3, assuming that the objects are spherical and 
that their dark matter distribution is given by the NFW 
density profile. The fit is done by minimizing \ 2 f° r the 
20 'data points' of <ji os and rei os together (the data points 
are independent because the number of particles in each ra- 
dial bin is very large, s ee the discussion in the Appendix of 
lLokas fc Mamonll2003l) . The data points were weighted by 
the assigned bootstrap errors although they do not account 
for the real variability of the data and therefore the quality 
of the fits in terms of \ 2 ls very bad. The best-fitting veloc- 
ity moments found for halo 1 are shown as dotted lines in 
Figures ISJEl 

In reality, when dealing with real galaxy data for clus- 
ters, the bins would include a few tens of objects instead 
of thousands, resulting in errors larger by at least an or- 
der of magnitude. To see how the sampling errors affect the 
determination of the dynamical parameters of clusters, we 
have also measured for each halo the velocity moments for 
a set of 40 randomly chosen particles per distance bin. This 
number of particles is chosen to be similar to the usual num- 
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Figure 10. Predicted profiles of <Ji os (R) and K\ oa (R) for dark 
matter halo with concentration c = 5.6 and different values of 
fi = —0.3, 0, 0.3, 0.6 and 0.9. The values of cr\ oa are scaled with 
circular velocity at 7*100 and distances are in units of rioo- These 
theoretical profiles are independent of the virial mass Mioo- 



ber of galaxies used in real clusters e.g. bv lLokas fc Mamonl 
( 2003 ). Errors were assi gned as described in the Appendix of 
lLokas fc Mamonl l)2003l) . We have then performed the same 
fitting procedure as described above. 

The best-fitting parameters for the haloes estimated 
with all dark matter particles are shown in the left column of 
Figure ITT1 as different symbols depending on the direction of 
observation with respect to the major axis: 0° (circles), 45° 
(triangles) and 90° (squares) . The crosses mark the 'real' val- 
ues of the parameters listed in Table Q In the right column 
of Figure im we show the best-fitting parameters estimated 
with 400 particles (40 particles per bin). 

Due to a rather time-consuming integration involved 
in the calculation of kurtosis, we restricted our analysis 



Table 2. Results of the fitting procedure 



Particles 


A log 


Afioo 


A log c 


A/3 


Alog((T r /cr 9 ) 


per bin 


mean 


Cr 


mean cr 


mean cr 


mean cr 


All 

10 


-0.03 
-0.07 


0.09 
0.10 


0.20 0.18 
0.08 0.24 


-0.78 1.04 
-0.20 0.48 


-0.12 0.12 
-0.04 0.11 



to the ten most massive haloes. The general conclusion is 
that when taking into account all three fitting parameters, 
for all haloes their best-fitting virial mass Mioo is reason- 
ably well recovered: the discrepancy between the best-fitting 
value and the real Mioo (measured using 3D information) is 
smaller than 62% (this large discrepancy is obtained for halo 
10 when observed at 45° with respect to the major axis). 
The next biggest error was obtained for halo 5 (37% when 
observed at 45° with respect to the major axis). 

In the case of the two remaining parameters, concen- 
tration c and anisotropy f3, all haloes show common trend 
in the discrepancies: /3 is underestimated while c overesti- 
mated, especially for the fitting including all dark matter 
particles. The discrepancies can be traced to the specific 
behaviour of the eri os and Ki os profiles. For example, in the 
case of halo 1 the three kurtosis profiles (see Figures 15171 
give similar /3 values, however the <ti os is much shallower 
for the observation angle of 45° than in the two remaining 
cases. As already mentioned, ki os is mainly sensitive to the 
velocity anisotropy. Since the line-of-sight velocity disper- 
sion profile can be made steeper in the centre either with a 
steeper density profile or with more radial orbits, and since 
the anisotropy parameter is almost the same in all three 
directions, we can expect the inferred concentration to be 
somewhat lower for the 45° direction, which is indeed the 
case. The situation is similar for halo 4. The kurtosis forces 
the anisotropy to be very tangential, so the very steep o"i oa 
profile at 0° (the leftmost middle panel of Figure [SJ requires 
a very large concentration of the density profile. 

We have calculated the means and standard deviations 
of the differences between the fitted parameters (logAfioo, 
log c, and f3 or log oy/o"e) and their values measured from the 
full 6D particle phase space. As expected from Figure lTTI we 
have found no significant dependence on the viewing angle. 
We show these quantities in Table [5] for the 30 cases stud- 
ied (3 for each of the 10 haloes) . We can see from this Table 
that, with 40 particles per bin, the virial mass is a little more 
underestimated, while the concentration parameter is less 
overestimated and the anisotropy less underestimated. This 
behaviour of c and f3 can be understood by recalling that 
in the case of using only 40 particles per bin the sampling 
errors of the moments are much larger. For very negative 
kurtosis small errors may enforce low f3 estimates, which 
have to made up by high c values (to reproduce velocity 
dispersion). Moreover, we find that the error bars listed in 
Tableware similar in both cases, which indicates that the 
physical variations due to substructure and to the different 
shapes of haloes dominate over the statistical noise or sam- 
pling errors. 

Table [3] shows the statistical significance of the biases 
on the parameters, using the Student's f-statistic, which as- 
sumes Gaussian distributions of the parameters (log Mioo, 
logc, (5 and ov/erg), and the binomial statistic testing the 
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Figure 11. Fitted (on projected data) values of virial mass Mioo, concentration parameter c and anisotropy j3 of the ten haloes for the 
three directions of observation with respect to the major axis of each halo: 0° (circles), 45° (triangles) and 90° (squares). The values 
measured on 3D data are shown with crosses. Left: Fitted parameters obtained with all particles. Right: Fitted parameters obtained 
with 40 particles per bin. 



distribution of the signs of these parameters. The two statis- 
tical tests lead to different conclusions about the bias of the 
parameters: the Student's t-statistic always leads to signifi- 
cant biases, while the binomial statistic indicates no signifi- 
cant bias in concentration parameter and velocity anisotropy 
when only 40 particles are used per radial bin. This differ- 
ence is caused by the skewed distributions of A log c, A/3 
and A log(ay/<jg). Although the Student's t-statistic has the 
advantage of being more sensitive to the outliers in the dis- 
tribution, it has the disadvantage of only being valid for 



Gaussian parent distributions, which is not the case here. In 
summary, in the case of 40 particles per radial bin, while 
the virial mass is biased towards lower values, we can- 
not conclude that the concentration parameter and velocity 
anisotropy are biased with the data we have. 

Given the range of Alogc seen in the right plots 
of Figure 1111 ( see al s o Tab le |5J, the discrepancies noted 
by lLokas fc Mamonl i2005tl between their concentration 
parameter for Coma an d t he lower values obta ined by 
iBiviano fc Girardil fcOO&t and lBiviano et al.1 (l2003f) in their 
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Table 3. Statistical significances in the biases of the fitting pro- 
cedure 



Particles 


log Afioo 


lo 


g c 


& 


log (ov/cr 9 ) 


per bin 


Pt P b 


Pt 


Pb 


Pt Pb 


Pt Pb 


All 


0.96 0.90 


1 


1 


1 1 


1 1 


10 


1 1 


0.96 


0.71 


0.96 0.82 


0.99 0.82 



Notes: Pt and Pt, are the probabilities for bias using the Student's 
t and binomial statistics, respectively. 



kinematical analysis of stacked clusters are reduced. As- 
suming that the departures in logc are independent of 
c and Mioo, we deduce that the difference in logc be- 
tween the measu rement for the Coma cluster (c = 9.4) by 
lLokas fc MjmrorJ and t he smaller concentration (c = 5.5) 
found bv lBiviano et alJ can be accounted for in 27% of our 
30 projected haloes. On the other hand, the larger difference 
in log c betwee n the Coma measurement and the concentra- 
tion found by iBiviano fc GirardU <|2003T) can be accounted 
for in only 10% of our projected haloes. Similarly, the differ- 
ence in logc between Coma a nd the value (c = 6) extrapo- 
lated from the values found bv lBullock et al.l J200ll) in their 
cosmological simulations, can be accounted for in 10 out of 
our 30 projected haloes. Moreover, there is a non-negligible 
scatter in the relation between h alo concentration and mass 
Jjing 1120001 : iBullock et al"ll200ll) . which reduces even more 
any discrepancy with high concentration found for Coma. 



5 DISCUSSION 

We studied the dynamical and kinematical properties of 
dark matter haloes obtained in cosmological A-body simu- 
lations. First, using all the 3D information available, we cal- 
culated their virial masses, radii, anisotropy parameters and 
estimated their density profiles. Next, we obtained projected 
velocity moments, standard observables used to estimate the 
dark matter content of virialized objects. We then fitted 
those 'data' with spherical models based on Jeans equations 
in order to reproduce the 'observed' velocity moments. 

Our approach was similar to the one of iTormen et alJ 
il997t) who used lowest order Jeans equation to model the 
velocity dispersion and find the masses of simulated haloes 
and compared them to the real masses of those haloes. 
In addition to velocity dispersion, we used the projected 
ku rtosis profiles in a simil ar way to that applied recently 
by lLokas fc Mamonl (1200311 to infer the properties of the 
Coma cluster. The use of kurtosis allows us to estimate 
the anisotropy of the velocity distribution. Besides the virial 
mass and anisotropy parameter, we fitted the concentration 
of the density profiles of the haloes. 

Our results emphasize the difficulties in the use of 
higher velocity moments to infer the properties of dark mat- 
ter haloes. The kurtosis seems to be very sensitive to the 
substructure and local matter flows. The discrepancies in 
the fitted properties of halo 4 can be traced to its peculiar 
mean radial velocity inside the virial radius which amounts 
to 0.26 fioo, while it is smaller than 0. If 100 for most of 
the remaining haloes (also in the tangential directions). We 
therefore confirm the necessity of using only elliptical galax- 
ies as tracers in the analysis of single clusters in order to 
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Figure 12. Line-of-sight velocity dispersion (top) and kurtosis 
(bottom) vs. projected distance to the centre measured at 0° to 
the principal axis of the ten most massive haloes in the simulation 
box. The thick line shows the results for halo 4, while the thin 
lines show the other nine haloes. 

minimize the effects of infall. Another source of problems 
lies in our very simple modelling of velocity anisotropy with 
constant f3, while this quantity really shows some radial de- 
pendence. Again, in the case of halo 4 this dependence is 
rather unexpected, departing from isotropy in the very cen- 
tre of the halo. 

In order to check to what extent the rather discrepant 
results for halo 4 are an exception or commonplace, we show, 
in Figure H2l the line-of-sight velocity dispersion and kurto- 
sis of the ten most massive haloes in the simulation box for 
a observer placed along the major axis of each halo. The 
thick line shows the results for halo 4, while for the rest of 
the haloes thin lines are used. We infer from Figure IT2l that 
the velocity moments of halo 4 follow an uncharacteristic 
pattern, probably caused by the unusually large negative 
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mean radial velocity of this halo (Table 0. We have also 
checked the pattern of the even velocity moments for the 
other two directions of observation used in Sections 3 and 4. 
In these cases, we found that the halo departing most from 
the general trend was halo 10, which has the highest mean 
velocity inside the virial radius with respect to Wioo (with a 
smaller ratio than halo 4). Therefore, local matter flows, as 
witnessed by non-zero mean velocity profiles, produce per- 
turbed line-of-sight velocity dispersion and/or kurtosis pro- 
files, which themselves can lead to an inaccurate estimate 
of the mass, concentration and/or velocity anisotropy of a 
cluster of galaxies. It is important to note, however, that 
concentration and anisotropy seem to be more affected than 
the mass estimation, which is quite robust. 

Figure Il2l also shows a general trend of the behaviour 
of the kurtosis of the simulated haloes. We notice that the 
kurtosis tends to be positive near the centre of the halo and 
negative at large distances. This means that the velocity dis- 
tribution is more peaked at the centre and more flattened 
outside than the purely Gaussian distribution in agreement 
with a recent finding of iKazantzidis. Magorrian fc Moore I 
(2003). But in spite of these departures from Maxwellian 
velocities at all radii, t he Jeans analysis used here and in 
lLokas fc Mamonl i2003T) produces fairly accurate measure- 
ments of the virial mass and concentration parameter (see 
Fig. [TTland Table 0. 

Although obvious interlopers to the haloes hav e been 
removed in a similar fashion as in lLokas fc Mamonl {2003, 
see also iKent fc Gurinlll982D . our results could depend on 
the neighborhood of the analyzed haloes. This would be the 
case if other haloes with similar velocity were present in 
the direction of observation, so that their particles would 
not be removed by the procedure described in Section 3. 
To estimate the plausibility of such a case, we have calcu- 
lated for every halo with mass up to one tenth of the mass 
of the most massive halo in the simulations (that makes 
a total of 40 haloes) the probability that choosing an ob- 
server at a random position around it, there would be one 
or more perturbing haloes with a mass greater than 25% of 
the mass of the studied halo, at a projected distance from 
the halo smaller than rioo and with a mean line-of-sight ve- 
locity within the interval [«i os — 3fioo,vi os + 3i>ioo]> where 
vios and vioo are the mean line-of-sight velocity and circular 
velocity at rioo of the analyzed halo. Choosing randomly 100 
observers for each halo, we found that the probability of not 
finding any perturbing neighbouring halo was higher than 
95% for 93% of the haloes, while for two of the analyzed 
haloes this probability was 88% and only one among the 40 
haloes had this probability smaller than 80% (77%). We can 
conclude from these numbers that it is quite unlikely that 
the results of the Jeans analysis are affected by neighbour- 
ing clusters. Moreover, the least isolated, i.e. the halo with 
the highest probability of having a perturbing halo around, 
among the top 10, halo 10, shows no specific bias in the 
parameter estimates (see Fig. II IB . 

Therefore, on one hand, the observed cosmic variance 
of the inner structure and internal kinematics of the mas- 
sive haloes in the cosmological simulations suggests that the 
typical properties of dark matter haloes are best obtained 
through the analysis of stacked observations as performed by 
ICarlberg. Yee fc Ellingsonl <ll997^lB"iviano fc Girardil (120031) 
and lBivian^^tHuTT2003F !But . on the other hand, in struc- 



tures (clusters of galaxies) with near zero mean velocity pro- 
files, this cosmic variance is much reduced. Therefore, it is 
well worth analyzing a single cluster with a large number of 
ve locities and a nea r zero mean velocity, such as was done 
bv lLokas fc Mamonl ll2003h for the Coma cluster. 
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